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Abstract 

Anisotropic  Cu-O  islands  form  striped  periodic  supergratings  upon  annealing  an 
oxygen-exposed  Cu(  110)  surface.  While  the  formation  of  these  striped  phases  has  been  attributed  to 
long-range  repulsive  forces  between  stripes,  we  show  using  the  one-dimensional  Ising  model  and 
Monte  Carlo  simulations  diat  short  range  adsorbate-adsorbate  interactions  are  capable  of  producing 
striped  phases.  In  these  systems  there  are  (at  least)  two  time  scales  of  importance;  the  time  required 
for  the  formation  of  Cu-O  strings  and  the  time  required  for  the  formation  of  stripes.  It  is 
demonstrated  that  the  latter  process  can  be  extremely  slow  at  room  temperature  and  that 
approximately  evenly  spaced  stripes  are  fonned  upon  annealing  by  employing  only  short  range 
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L  bitrodncik» 

Striped  phases  are  observed  in  nature  for  diverse  systems  from  sorted  stone  stripes  to  clouds  [  1  ] . 
Kem  et  al  observed  that  when  oxygen  is  adsorbed  on  Cu(llO)  and  annealed,  anisotropic  Cu-O 
islan^ts  form  Striped  periodic  supeigratings  [2].  The  observed  spacings  between  Cu-O  islands 
depended  upon  the  coverage  and  temperature,  and  were  studied  using  scanning  tunneling 
microscopy  (STM)  and  He  diffraction  [2,3].  The  unannealed  0/Cu(  110)  surface  has  been  the  subject 
of  a  number  of  earlier  STM  smdies,  and  displays  disordered  “strings”  of  adsorbates  along  the  <00 1  > 
direction  which  are  one  atom  wide  in  the  <lT0>  direction  [4-6],  or  islands  of  condensed  strings 
[5,7].  These  “added  row”  strings  are  found  on  amimber  of  O  and  S  covered  fcc(l  10)  surfaces  [8,9]. 
Kero  et  al  found  that  annealing  causes  several  strings  to  condense  together  forming  “stripes”  of 
quite  uniform  width  and  spacing  in  the  <00 1>  direction.  These  interesting  results  have  been 
interpreted  as  being  the  result  of  long  range  forces  due  to  surface  stresses  [2,10].  However,  such 
forces  alone  are  unable  to  account  for  the  narrow  distributions  of  stripe  widths  and  spacings  found 
experimentally  [10-12]. 

In  this  paper  we  use  a  one-dimensional  approximation  and  two-dimensional  simulations  to 
explore  the  roles  of  short  range  forces  and  interface  entropy  (or  the  so-called  “entropic  repulsion" 
[I3.I4|)  on  the  formation  of  superstructures.  In  the  limit  of  a  highly  asymmetric  interaction,  a 
two-dimensional  lattice  can  be  approximated  by  a  one-dimensional  lattice.  Given  the  strong  Cu-O 
interaction  in  the  <00 1  >  direction,  we  use  the  one-dimensional  Ising  model  to  approximate  the 
0/Cu(l  10)  system  in  the  <lT0>  direction.  For  this  simple  and  well  known  model  we  show  that 
evenly  spaced  islands  are  the  most  likely  configuration.  Additionally,  we  have  performed  a  Monte 
Carlo  study  of  the  two-dimensional  .system.  The  experimental  coverage  dependence  on  stripe 
spacing  for  CVCu(  1 10)  are  reproduced  in  both  cases. 

IL  One-Dimcnsioiial  Ising  Model 

For  simplicity,  the  one-dimensional  Ising  model  is  discussed  in  terms  of  the  micro-canonical 
ensemble;  that  is.  there  are  N  atoms  on  M  sites  with  N 1 1  nearest  neighbors  [15].  The  energy,  which  is 
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set  equal  to  tbe  ejqpectatiou  value  of  llie  eneiigy  for  die  correspo&dmg  canonical  system,  can  be 
expressed  as 

AN-  NJ2)  (1) 

wlieie  Nqi  is  die  number  of  nearest  neighbor  pairs  widi  one  site  fiiU  and  one  site  empty  and  J  is  the 
interaction  eneigy.  Thus,  the  enetgy  is  only  dependent  on  the  number  of  islands.  The  most  likely 
distribution  of  atoms  amtuig  the  islands  is  the  one  with  the  highest  degeneracy.  The  number  of  ways 
of  arranging  the  atoms  into  the  islands  so  that  there  are  ni  in  the  first  island,  n2  in  the  second,  etc.  is 
given  by 

- 

^  "  n,!n2!n3!..7i^^!  (2) 


(where  N's  NoNoi).  which  is  maximum  for  ni »  02  ^  03  s ...  s  dn.  A  similar  argument  holds  for  the 
empty  sites.  Thus,  the  most  likely  distribution  has  islands  of  equal  thickness  and  equal  spacing. 

Fig.  1  shows  the  most  probable  island  spacing  as  a  function  of  temperature  and  coverage  for  the 
micro-canonical  one-dimensional  Ising  model.  The  functional  forms  of  these  curves  are  in  good 
agreement  with  the  experimental  observations  for  annealed  0/Cu(1 10).  However,  the  expected 
fluctuations,  even  for  the  micro-canonical  ensemble,  are  large.  Below  it  is  demonstrated  that 
asymmetric  short  range  interactions  on  a  two-dimensional  lattice  can  give  evenly  spaced  stripes 
with  high  probability. 

in.  IWo-Diniensional  Monte  Carlo  Simulations 

The  Monte  Carlo  simulations  performed  here.  (16].  use  similar  interaction  energies  to  those 
reported  earlier  (6]  with  the  important  addition  of  a  next-nearest  neighbor  attraction  in  the  <lT0> 
direction  as  suggested  by  VTitterlin  et  al  (7  J.  Additional  interactions  include  stability  for  strings  of  at 
least  four  adatoms.  This  is  consistent  with  the  STM  data  in  which  strings  of  several  adatoms  were 
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dMerved  to  move  as  umts  across  die  surfiikce  [4-6].  The  interacticMiaieipestfiat  we  use  are  shown  in 
Fig.  2  and  are  described  below. 

Aldiough  O  and  Cu  are  both  mobile,  it  is  sufficient  to  consider  one  mobile  species  in  the  Monte 
Cario  gmiilarinns.  These  are  referred  to  below  as  nuclei,  and  they  consist  of  an  oxygen  adatom 
bonded  to  a  single  Cu  adatom  with  die  bond  between  the  two  unbreakable.  The  nuclei  are  adsorbed 
on  a  fcc[110]  lattice.  The  number  of  nuclei  used  was  determined  by  the  experimental  coverage  of 
oxygen  (0)  and  was  conserved  for  the  duration  of  the  simulation  (canonical  ensemble).  The 
diffusion  of  the  nuclei  was  modeled  as  isotropic.  The  short  range  interactions  used  in  this  model 
were  initially  approximated  by  the  behavior  of  this  system  as  observed  by  STM  [4-6],  and  where 
originally  used  to  describe  the  structure  of  the  unannealed  surface  [6].  First,  the  <00 1> 
nearest-neighbor  interaction.  Fig.  2a.  was  chosen  to  be  attractive  as  is  implied  by  the  formation  of 
strings.  Next,  the  <(X)1>  diagonal  next-nearest  neighbor  interaction.  Fig.  2b,  is  included  to  account 
for  die  experimental  observation  that  kinked  strings  are  also  commonly  observed  [6].  The  attractive 
interaction  between  strings  along  the  <1I0>  direction.  Fig.  2c.  is  included  to  reproduce  the 
observation  of  coalescence  of  strings  into  stripes  and  its  magnitude  is  on  the  order  of  the  interaction 
energy  measured  by  Besenbacher  er  ai  [17|.  The  only  repulsive  interaction  is  the  <110> 
nearest-neighbor.  Fig.  2d.  This  interaction  is  invoked  to  give  the  void  or  missing  row  observed 
experimentally,  and  the  preferential  growth  of  Islands/stripes  in  the  <(X)  1  >  direction.  The  additional 
four  nuclei  chain  interaction.  Hgs.  2e  and  2f.  accounts  for  the  observation  that  chains  of  four  or  more 
nuclei  are  more  stable  than  chains  of  fewer  nuclei  and  move  across  the  surface  as  a  linked  unit  [5]. 

Starting  with  randomly  distributed  nuclei  at  640K.  we  find  that  strings  of  nuclei  rapidly  form. 
Upon  cooling  from  640K  to  293K  fairly  evenly  spaced  stripes  of  nuclei  form.  Typical  Monte  Carlo 
simulation  results  are  shown  in  Fig.  3.  The  numbers  and  spacings  of  these  stripes  depend  upon  the 
coverage,  but  are  the  same  for  multiple  runs  at  a  given  coverage.  We  find  that  from  ca.  8x10^  up  to 
ecu  10^  moves  per  nucleus,  we  obtain  identical  results  for  lattices  of 30x30 Cu  atoms.  We  find  similar 
results  for  lattices  of  30x30. 60x60.  and  128x128  Cu  atoms,  as  well  as  for  asymmetric  lattices,  all 
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with  periodic  boundary  conditions.  The  structures  formed  were  also  independent  of  the  initial 
configuration.  The  same  striped  phases  resulted  either  from  nuclei  initially  distributed  randomly  or 
when  one  large  complete  island  (stripe)  initially  contained  all  nuclei. 

Fig.  4  shows  the  simulated  stripe  spacing  (D)  vs.  coverage  using  only  the  short-range 
interactions  given  above.  This  result  is  in  qualitative  agreement  with  the  experimental  OATufl  10) 
data  —  Fig.  2a  of  Ref.  [2].  The  small  fluctuations  in  our  data,  obtained  by  several  repeated 
Qifniilatinns  at  the  Same  coverage,  are  q)proximately  the  size  of  the  points  in  the  Fig.  4.  Note  that  the 
vertical  scale  of  the  data  in  Fig.  4  varies  dramatically  with  the  interaction  energies  chosen. 

In  our  simulations,  by  randomly  adsorbing  Cu-O  nuclei  at  294K  followed  by  the  same  number 
of  steps  as  for  high  temperature/aiuiealed  adsorption,  we  reproduce  the  previous  STM  data  [4-6]  and 
simulations  [6]  which  resulted  in  disordered  strings.  These  starting  conditions  resulted  in  some 
clustering  of  strings  as  seen  by  STM  in  Ref.  [7J.  However,  starting  with  one  perfect  stripe  several 
nuclei  wide,  we  find  that  at  294K  in  10*^  attempted  moves  per  nuclei  that  this  island  broke  up  into  the 
same  number  of  stripes  found  for  the  annealed  surface  at  this  coverage.  We  take  this  as  evidence 
pointing  toward  the  striped  phases  being  equilibrium  structures  in  our  simulations. 

iV.  Condusions 

By  con.sidering  the  strong  adatom  attraction  in  the  <0()1>  direction  as  the  determining 
interaction  in  string  formation,  and  the  weaker  adatom  attraction  in  the  <lT0>  direction  as 
determining  the  string-string  interactions,  we  modeled  the  O/Cu(110)  system  using  a 
one-dimensional  Ising  model.  Even  in  this  simple  equilibrium  model  equally  spaced  bands  were 
obtaiited  with  spacings  and  widths  i*s.  coverage  and  temperature,  in  qualitative  agreement  with  Fig.  2 
of  Ref.  12). 

The  es.sential  qualitative  features  observed  in  the  STM  measurements  of  [2-6]  on  the  0/Cu(  1 1 0) 
surface  at  room  temperature  and  upon  annealing  have  been  reproduced  using  simulations  which 
consider  only  short  range  interactions  including  a  weak  adatom  attraction  in  the  <1 10>  direction. 
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This  is  in  cmiinst  to  the  suggestion  that  loog  range  forces  need  be  invoked  in  order  to  explain  the 
antwaWiiH  data  [2,18].  It  has  been  demonstrated  diat  long  range  ordering  can  be  obtained  with  short 
range  interactions.  Tliis  behavior  is  evident  in  the  one-dimensional  Ising  model,  albeit  with  laige 
fluctuations.  In  two-dimensional  Monte  Carlo  simulations  the  behavior  exists  with  much  smaller 
fluctuations. 

The  formation  of  the  equally  spaced  stripes  on  die  (VCu(l  10)  surface  is  at  least  partially  due  to 
that  particular  configuration  simply  having  the  highest  degeneracy  (from  the  Ising  model).  The 
interface  entropy  leads  to  distributions  of  stripes  with  approximately  equal  widths  and 
approximately  equal  spacings.  We  believe  that  the  interface  entropy  and  short  range  forces  play  an 
important  role  in  the  formadon  of  the  experimentally  observed  striped  supergradng  phases  for 
CVCu(l  10).  As  with  stress,  short  range  forces  and  entropy  alone  might  be  insufficient  to  describe  the 
narrow  experimental  distribudons  of  stripe  widths  and  spacings.  This  would  be  analogous  to 
step-step  interacdons  in  which  Bartelt  tt  at.  have  shown  that  the  narrow  distribudon  of  terrace 
widths  on  surfaces  are  due  to  both  entropy  and  stress  terms  [  1 3, 1 4].  We  conclude  that  a  combination 
of  enuopy  and  stress  should  lead  to  the  observed  narrow  stripe  distribudons  for  0/Cu(  1 1 0).  Future 
direcdoru;  include  modeling  the  temperature  dependence  and  the  kinedcs  of  the  formadon  of  this 
system  as  well  as  other  possible  causes  of  long  range  ordering  [  19]. 
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nCURE  CAPTIONS 

1 .  Stripe  q)acing  (D)  as  a  function  of  a)  coverage  (at  J/Ts  -5)  and  b)  temperature  at  half  monolayer 
coverage  for  the  1-D  Ising  model  with  M=60  sites. 

2.  Interactions  used  in  the  Monte  Carlo  simulations.  The  second  layer  Cu  atoms  of  the  Cu(l  10) 
surface  t^pear  as  die  large  white  circles,  the  top  layer  Cu  atoms  appear  as  dark  gray  circles,  the  O 
atoms  appear  as  small  white  circles,  and  the  reference  CuO  nuclei  appear  in  solid  black.  The 
black  heavy  ringed  top  layer  Cu  atoms  highlight  the  CuO  nuclei  interacting  with  the  reference 
nuclei:  a  solid  white  with  black  ring  denotes  one  possibility,  while  gray  with  black  ring  denotes 
all  equivalent  possibilities.  The  stability  gained  by  the  reference  nuclei  is  in  each  case:  a)  -109 
meV,  b)  -109  meV,  c)  -26  meV,  and  d)  +1690  meV.  In  addition,  the  stability  gained  if  there  are 
four  nuclei  in  a  string  is:  e)  -77  meV  and  f)  -154  meV, 

3.  Typical  examples  of  simulation  results  for  coverages:  a)  0. 1 6  ML,  b)  0.22  ML,  c)  0.26  ML  and 
d)  0.37  ML.  Periodic  boundary  conditions  are  used. 

4.  Stripe  spacing  (D)  vs.  O  atom  coverage  on  the  Cu(  1 1 0)  surface  after  cooling  from  640K  to  294K 
deduced  from  Monte  Carlo  simulations  allowing  ca.  1 0^  moves  per  adatom.  The  displayed  data 
point  size  approximately  represents  the  error. 
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